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ABSTRACT 

Observations have demonstrated the presence of substantial amounts of interstellar 
dust in elliptical galaxies, most of which is believed to be distributed diffusely over 
the galaxy. Absorption by dust grains has a major impact on the transfer of stellar 
radiation, and affects the projection of each physical (and kinematic) quantity. In a 
previous paper, we have investigated the effects of a diffusely distributed dust compo- 
nent on the observed kinematics of spherical galaxies. In this paper we investigate the 
effect of not taking dust into account in dynamical modelling procedures. We use a set 
of semi-analytical dusty galaxy models to create synthetic LOSVD data sets, which 
we then model as if no dust were present. 

We find some considerable differences between the best-fitting models and the 
input models, and we find that these differences are dependent on the orbital structure 
of the input galaxy. For radial and isotropic models on the one hand, we find that the 
dynamical mass of the models decreases nearly linearly with optical depth, with an 
amount of 5 per cent per optical depth unit, whereas their orbital structure is hardly 
affected. For tangential models on the other hand, the dynamical mass decrease is 
smaller, but their orbital structure is affected : their distribution functions appears less 
tangentially anisotropic. For all models the mass-to-light ratio will be underestimated, 
by a factor of around 20 per cent per optical depth unit. 

We discuss these results in the light of the limited effects of dust extinction on 
the LOSVDs, as obtained in paper I, and conclude that the determination of the 
dynamical mass and the kinematic structure of galaxies is not only determined by the 
observed kinematic quantities, but is also critically dependent on the potential and 
hence the observed light profile. We argue that dust, even in rather modest amounts, 
should therefore be taken into account in kinematic modelling procedures. 

Key words: dust, extinction - galaxies : elliptical and lenticular, cD - galaxies : ISM 
- galaxies : kinematics and dynamics 



1 INTRODUCTION 

It has become well-established that early-type galaxies con- 
tain a considerable amount of interstellar dust. In the opti- 
cal dust is detected in ellipticals by its obscuration effects 
on the light distribution, when it is present in the form of 
dust lanes and patches (e.g. Ebneter & Balick 1985, Veron- 
Cetty & Veron 1988, van Dokkum & Franx 1995). In emis- 
sion, dust is detected by the IRAS satellite in the 60 and 
100 (im wavebands (Knapp et al. 1989, Roberts et al. 1991). 
In a comparative analysis Goudfrooij & de Jong (1995, here- 
after GdJ95) show that the dust masses derived from FIR 
data are about a factor of ten higher than those calculated 
from optical data. Using a more detailed dust mass estimator 
which includes a temperature distribution for the dust grains 



(Kwan & Xie 1992), Merluzzi (1998) shows that GdJ95 still 
underestimated the FIR dust masses with a factor up to 
six. Also submillimeter observations (Fich & Hodge 1991, 
1993, Wiklind & Henkel 1995) and ISO data beyond 100 fira 
(Haas 1998), that may be able to detect the very cold dust 
for which IRAS is insensitive (T < 25 K), suggest that the 
dust masses could be up to an order of magnitude higher 
than observed from FIR observations alone. 

This difference between the absorption and emission 
dust masses in elliptical galaxies is called the dust mass dis- 
crepancy. It cannot be solved by a more critical reconsidera- 
tion of the IRAS data (Bregman et al. 1998), by corrections 
for the optical absorption in the dust lanes (Merluzzi 1998) 
or by taking into account the dust recently ejected from 
evolved stars (Tsai & Mathews 1996). The interstellar dust 
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medium in ellipticals hence has to be composed of (at least) 
two components: a lesser massive one which is optically visi- 
ble in the form of dust lanes and patches, and a more massive 
one which is distributed over the galaxy and is hard to de- 
tect optically. The spatial distribution of this component is 
still unclear. GdJ95 suggest that it is distributed over the in- 
ner parts of the galaxy, which supports the evaporation flow 
scenario. In this picture, most of the gas and dust has an 
external origin: clouds of interstellar matter have been ac- 
creted during interactions/merging with a gas-rich galaxy, 
and they gradually evaporate in the hot X-ray emitting gas 
(Sparks, Macchetto & Golombek 1989, de Jong et al. 1990, 
Forbes 1991, GdJ95). Tsai & Mathews (1996) and Wise & 
Silva (1997) on the other hand suggest that dust is not con- 
fined to the inner few kpc only, but also extends to larger 
radii. Such a dust distribution supports a scenario where 
gas and dust have an internal origin: it can be associated 
with red giant winds or formed in connection with star for- 
mation in cooling flows (Fabian, Nulsen & Canizares 1991; 
Knapp, Gunn & Wynn- Williams 1992; Hansen, j0rgensen 
& N0rgaard-Nielsen 1995). However, none of these mech- 
anisms seems to be able to explain the FIR observations 
satisfyingly, and a combination of both mechanisms may be 
at work (Merluzzi 1998). As in the first paper of this series 
(Baes & Dejonghe 2000, hereafter paper I), we will assume 
that dust is distributed smoothly over the entire galaxy. 

Absorption by dust grains has a major impact on the 
transfer of stellar radiation through the interstellar medium. 
It therefore affects the projection, i.e. the integration along 
the line-of-sight (hereafter LOS) of the intrinsic three- 
dimensional light distribution. A number of studies have 
investigated the effects of a diffuse dust distribution upon 
the photometry of ellipticals (Witt, Thronson & Capuano 
1992, GdJ95, Silva & Wise 1996, Wise & Silva 1996; here- 
after WS96). But dust does not only affect the projection of 
the light distribution, it affects the projected kinematics too. 
This is of particular importance in stellar dynamics, where 
the ultimate purpose is the determination of the phase space 
distribution function F(r,v), hereafter DF, describing the 
entire dynamical structure. This DF is usually constructed 
by fits to the LOSVDs or their moments, and hence it is 
obviously important to examine to which extent these are 
affected by dust obscuration. 

In paper I we constructed a set of semi-analytical mod- 
els in order to investigate the effects of dust absorption on 
the light profile, the projected velocity dispersion profile and 
the LOSVDs. Concerning the photometry, we find that dif- 
fuse dust has a strong impact : a global attenuation, strong 
extinction in the central regions, the formation of radial 
color gradients, and an increasing apparent core radius. In 
spite of the simplified character of our analysis (we don't 
include scattering effects in our models), these results are in 
qualitative correspondence with the conclusions of the more 
detailed photometric studies mentioned before. The effect of 
dust on the projected dispersion profile or the LOSVDs is 
of a totally different nature. The effects are a redistribution 
along the LOS: dust makes the contribution of the nearer 
parts more important, relative to the more distant parts. 
Therefore it is no big surprise that, in a spherically sym- 
metric galaxy, the LOSVDs are quite insensitive to modest 
amounts of dust. For example, for a modest optical depth 
(r = 2) the effect on the projected velocity dispersion is 



around 2 percent in the central regions. However, for higher 
optical depths, these effects do become considerable, which 
can be important on a local scale: e.g. asymmetries in the 
projected kinematics may be the result of a large extinc- 
tion in a compact region, such as in the dust lane elliptical 
NGC 5266 (Mollenhof & Marenbach 1986). More details can 
be found in paper I. 

The fact that both the photometry and the observed 
kinematics are affected by dust obscuration will have conse- 
quences on the dynamical modelling of galaxies. Now that 
we understand the way dust absorption affects the projected 
kinematics, it is obvious that dust needs to be accounted for 
in deprojection procedures, or more general, in kinematic 
modelling procedures. It is the scope of this paper to de- 
termine the effects of neglecting dust in the modelling of 
kinematic data. In other words, if one tries to model dust- 
affected kinematic profiles without taking dust into account, 
as is usually done, to what extent are erroneous conclusions 
drawn about the kinematic structure of the galaxy ? In order 
to answer this question we will create synthetic dust-affected 
data sets, and model them as if no dust were present. 

In section 2 and section 3 we describe the construction 
of the data sets and the modelling procedure. The results 
of the modelling are presented in section 4, and in section 5 
we investigate whether these depend critically on the choice 
of the dust geometry. A discussion of the results is given in 
section 6. 



2 THE DATA SETS 
2.1 The model 

To create our data sets, we construct a set of two-component 
galaxy models, consisting of a stellar and a dust component. 
Both components are spherically symmetric. 

For the stellar component, we use a Plummer model 
(Dejonghe 1987), which is described by the potential-density 
pair 



i/j(r) 



GM 



, , 3 Mo 



1 + 



(la) 



(lb) 



with c = 5 kpc the so-called core radius and Mo = 5 x 
1O 1O M the total mass. Furthermore we assume a con- 
stant mass-to- light ratio T(r) = 4 Y©, such that the three- 
dimensional light distribution is given by £(r) = jy^ p( r )- 
With the technique described by Dejonghe (1986) 
one finds a family of two-integral DFs F q (E,L) that 
self-consistently generate the Plummer potential-density 
pair (p"a|b). This family has two interesting properties, which 
justify its choice as a generic model for the class of elliptical 
galaxies (although its projected density profile does not fit 
real elliptical galaxies). First, the models are completely an- 
alytical, i.e. the DF and all the (projected) moments can be 
calculated analytically. And second, the family depends on 
one single parameter q which can be varied continuously in 
order to obtain tangential (q < 0), isotropic (q — 0) or radial 
(0 < q < 2) models. We consider a set of different orbital 
structures, characterized by the parameters q — —6, —2, 
and 1. In the sequel of this paper we will refer to them 
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Table 1. Some parameters of the dusty Plummer models as a 
function of the optical depth t. The second column gives the 
observed core radius, the third and the fourth column give the 
observed luminosity and the total extinction, and column six and 
seven list the projected light density and projected dispersion for 
the central LOS. This dispersion is tabulated for the isotropic 
model, one finds the central dispersion of the very tangential, 
tangential and radial models by multiplying this with the factors 
%/2/2, V3/2 and ^6/5 respectively. 



T 


^obs 


^obs 


A 


<?P,0 


a P,0 




(kpc) 


(10 9 L ) 


(mag) 


(L /pc 2 ) 


(km/s) 


0.0 


5.00 


12.50 


0.00 


159.2 


178.0 


0.5 


5.33 


11.08 


0.13 


124.7 


177.8 


1.0 


5.65 


9.91 


0.25 


99.0 


177.6 


1.5 


6.02 


8.93 


0.37 


79.5 


177.1 


2.0 


6.35 


8.11 


0.47 


64.6 


176.5 


2.5 


6.74 


7.41 


0.57 


53.1 


175.7 


.3.0 


7.15 


6.81 


0.66 


44.2 


174.9 



as very tangential, tangential, isotropic and radial models 
respectively. 

As in paper I, we only incorporate the effects of dust ab- 
sorption and neglect scattering effects. Then, the dust com- 
ponent is completely determined by the opacity function 
n(r), for which we use a modified Hubble profile 



re(r) 



t I r 
2^' 1 + - 



(2) 



The normalization is such that r equals the total optical 
depth, denned as the projection of the opacity along the 
entire central LOS^j 

r r+co 

t= n(r) ds = 2 / ft(r) dr. (3) 

./central LOS J 

As WS96 and paper I, we take the same core radius for 
the dust and the stars; the choice of the opacity function is 
critically investigated in chapter 5. We only consider optical 
depths ranging from r = to r = 3 in our calculations, as 
high values for the optical depth associated with a diffuse 
dust component seem not to be in accordance with photo- 
metric studies (GdJ95, WS96). 



2.2 The projections 

For each model we create a set of so-called dusty projected 
kinematic data. A dusty projected quantity fj, p (x,v) differs 
from a normal projected quantity as it is a weighted integral 
along the LOS a; of a three-dimensional spherically symmet- 
rical quantity n(r,v). We assume that the galaxy is located 
at a distance which is significantly larger than its size. Then 
the errors made by assuming parallel projection, which are 
of the order (c/D) 2 with D the distance to the galaxy, are 
negligible (see Paper I, section 2). For a dusty galaxy with 
opacity function ft(r) the appropriate formula reads 



H P {x,v) = 2 



K(x, r) 



fx(r,v)rdr 
\/r 2 — x 2 



(4a) 



* Throughout paper I and this paper we use this definition of T, 
whereas e.g. WS96 defined r as the integral of the opacity from 
the centre of the galaxy to the edge, half our value. 



where K[x, r) is a weight function defined as 
x n(r)rdr 



K.(x, r) — exp 



cosh 



re(r') r' dr' 
Vr' 2 — x 2 



(4b) 



Details can be found in section 2 of paper I. We substituted 
the analytical expressions for the moments of the DF, to- 
gether with the opacity function (jij), in expression ([la]) to 
obtain dusty projected profiles, such as the projected light 
density l v (x) and the projected velocity dispersion a p (x). 
Since in particular the projected light density will depend 
on the optical depth, we will have to be careful in our ter- 
minology when calculating quantities such as the luminosity 
and the core radius. We need to discriminate between true 
and observed quantities, which are respectively derived from 
the spatial and the projected distribution. For example, the 
true luminosity of the galaxy is calculated by integrating the 
light density over space 



L — 4-k 



£(r) r dr, 



(•>) 



which is of course independent of the optical depth and 
equals Mo/T — 1.25 x 10 10 Lq. The observed luminosity 
L bs is calculated by integrating the projected light density 
over the plane of the sky, 



L h 



2tt 



£ p (x) x dx, 



(6) 



and is a function of r. Analogously, all our Plummer models 
have a (true) core radius c = 5 kpc, whereas their observed 
core rad ius c b s will depend on the optical depth (see sec- 
tion 3.1). In table [j] we tabulate some of the parameters of 
our dusty Plummer models. 



2.3 The data sets 

Each data set consists of the projected light density data 
£ p (x) and projected dispersion data a p (x), which are taken 
at x — Okpc, 0.5 kpc, . . . ,10 kpc. However, since dispersion 
profiles depend on both the orbital structure and the mass 
distribution, they do not sufficiently constrain the dynami- 
cal structure of the galaxy. This mass-anisotropy degeneracy 
can be broken by including the higher order Gauss-Hermite 
moments in the fitting routine (van der Marel & Franx 1993, 
Gerhard 1993). Since we work with simulated data, we are 
able to include LOSVD data points directly in the modelling 
procedure. We assume that these data are noiseless, such 
that our data set is perfect, i.e. it contains all the kinematic 
information that could be available from perfect observa- 
tions. Altogether each data set consists of 252 data points. 



3 THE MODELLING PROCEDURE 

3.1 Determination of the potential 

The first step in the modelling of kinematic data is the de- 
termination of the potential. Although we have a perfect 
data set at our disposal, we are not able to constrain the 
potential completely without any assumptions, not even in 
the case of spherical symmetry (Dejonghe & Merritt 1992). 
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Figure 1. The difference A/i(x) in surface brightness between 
our dusty projected density profiles and the best fitting Plummer 
light profiles, for the modest optical depths r = 1 and r = 2, 
and also for the higher values r = 5 and r = 10. It is clear that 
for modest optical depths a Plummer profile provides a satisfying 
fit to the dust-affected profiles, such that we can safely assume a 
Plummer potential in our modelling procedure. The dots in the 
figure indicate the observed core radii c oba of the Plummer light 
profiles, where by construction ^ p (c obs ) = ^p'(c b s ). 



In this first approach we assume that it makes sense to con- 
sider a constant mass-to-light ratio, such that the spatial 
dependence of the potential can be derived from the pro- 
jected light density £ p (x). It can be calculated numerically 
for our set of models, but we prefer to work with a poten- 
tial whose functional form is explicitly known. We find that, 
for the modest optical depths we are considering, the dust- 
affected light profiles of our models can still be satisfyingly 
described by Plummer light profiles 



tl\x) 



1 + 



(7) 



with c bs dependent on the optical depth. In figure [l] we plot 
the difference 



Afi(x) = -2.5 log \l p (x)/lf{a 



(8) 



in surface brightness between the dusty projected light den- 
sity profiles l p (x) and the best fitting Plummer light profiles 
Even at very large projected radii we find that A/j, 
never becomes larger than 0.05 for r = 2, such that we can 
say that the Plummer character of the galaxy is preserved 
for modest values of r. We will therefore assume a Plummer 
potential 



GM 
ip(r) = 

Cobs 



l + ■ 



(9) 



for our models, with Cobs determined from the best fit to 
the £ p (x) data, and the mass M still featuring as a free 
parameter. The values for c b s are tabulated in the second 
column of table 1. 



& Merritt 1992). It is always possible to write this DF as a 
infinite sum of components 



F(E, I?) = ^c,F\E,L 2 



(10) 



where d are the coefficients and the components F'(E,L 2 ) 
form a complete set of simple dynamical models. For any 
observed kinematic data point fj, n (x,v) the same expansion 
is valid, 



= E 



(11) 



since these moments depend linearly on the DF. Practically 
one can only consider a finite number N of components, 



F(E, L 2 ) w ^ c t F'(E, L 2 ). 



(12) 



The best fitting coefficients can then be found by minimizing 
a x 2 -like variable 



X 



E 



EN i 



(13) 



where the sum contains all data points, and w n is the weight 
accorded to the nth data point. Since we assume to obtain 
perfect, noiseless synthetic data, we can use these weights to 
(arbitrarily) set the relative importance of each data point in 
the global x 2 . For the projected density points we take w n — 
1, for the projected dispersion w n — | and for the LOSVD 
data the weight varies from w n = | at the centre to w„ = 
in the outer parts. This x 2 is quadratic in the coefficients 
and has to be minimized under the linear constraint that the 
DF has to be positive over a grid (EjjL 2 ,) in phase space, 



E c ^ 1 



E 3 ,L 2 )>0 



for all j and k, 



(14) 



which amounts to a typical Quadratic Programming prob- 
lem. For details we refer to Dejonghe (1989). 

We choose our components from a library of Fricke mod- 
els. These are simple dynamical models that are defined by 
the augmented densityn, 



(15) 



where s* is a scale factor and a and 6 are real numbers that 
satisfy the condition a — 2b > 3 to keep the total mass finite. 
The anisotropy /3(r), generally defined as 



(16) 



will be constant for these models, since the r-dependence is 
a power law (Dejonghe 1986, section 1.5.1). We find immedi- 
ately /3(r) = —b, hence the condition b > — 1 is required. The 



3.2 Determination of the DF 

With a fixed potential there exists one and only one two- 
integral DF F(E, L 2 ) that fits the kinematic data (Dejonghe 



* The augmented density is the density written as a function of 
the radius r and the potential i/j. It is a fundamental quantity in a 
technique to construct 2-integral distribution functions (Dejonghe 
1986). 
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Figure 2. Plot showing \ 2 as a function of the mass M, for a 
tangential input model with r = 1. The x 2 values have no abso- 
lute meaning. The dynamical mass is determined at the minimum 
of the plot, being 4.85 X 10 10 Mq. The plot is shown for models 
containing 10, 15, 20 and 25 components. Clearly 20 components 
are sufficient to determine the dynamical mass accurately. 



two-integral DF corresponding to (|l^) is a simple power-law 
of E and L 2 



F(E, L 2 ) oc E a - b -iL 2b . 



(17) 



For tangentially anisotropic components the central density 
vanishes, for isotropic ones it is finite and non-zero and ra- 
dially anisotropic components have a central density cusp. 
The advantage of this family lies within the fact that most 
of the kinematics can be calculated analytically for positive, 
integer values of b (De Rijcke & Dejonghe 1998). 



3.3 Determination of the dynamical mass 

The only unknown in our model now is the mass M, which 
still acts as a free parameter. For its determination we run 
our models for a number of possible values, and determine 
the best fitting DF and the according \ 2 parameter for each 
value. The best fitting mass of the model is then determined 
as the minimum in x 2 (M). 



3.4 Practical application 

Practically we construct models for a sufficiently large set 
of masses, with a step AM = 5 x 10 8 Mq, 1 per cent of the 
input mass. For each mass we construct a DF using a li- 
brary of 30 isotropic and tangential Fricke components. For 
computational reasons we do not include any radial com- 
ponents; radial models can be constructed by linear com- 
binations of isotropic and tangential components, in which 
the latter have a negative weight (allowed as long as the 
DF remains positive over phase space). We typically used 
20 components: adding more components does not signifi- 
cantly affect the results (neither the DF nor the mass) any- 
more. This is illustrated in figure ^, where we plot the x 2 
values of the various fits in function of the number of com- 
ponents and the value of the total mass. 

A check on the fitting procedure is done by using r = 
data, i.e. data from non-dusty Plummer models. In order to 



make the test robust, we use the same strategy as De Rijcke 
& Dejonghe (1998): for the Plummer models that can be 
fit exactly in terms of Fricke components (q = —6, —2 and 
0), we remove these components from the library. Still, both 
the DF and the mass of the input models can be successfully 
reproduced. 

In figure ^ we plot some results of the fitting procedure 
for the four r = 2 models. Shown are the projected light den- 
sity £ p (x), the projected dispersion cr v (x) and the LOSVDs 
at x = kpc and x — 10 kpc, the innermost and outermost 
LOSVD in our dataset. The quality of the fit cannot be 
deduced from the x 2 values, since these have no statistical 
meaning. At bottom row in figure ^ we show the projected 
fourth moment (, p {x) of the DF, which is not included in the 
modelling procedure. It can be used to check the quality of 
the fit, which is very satisfactory in all cases. 



4 RESULTS 

In this section we describe the results of our fitting procedure 
and the kinematic properties of the models. For the sake of 
clarity we first explicitly define some terms. 

The input models are the models that are described in 
section 2.1, i.e. Plummer galaxy models containing a dust 



component, whereas the fitted models are the models that 
come out of the modelling procedure, and which, by con- 
struction, contain no dust. Since each couple of parameters 
(q, t) corresponds to one input model, and hence one data 
set and one fitted model, we will call the combination of 
input and fitted models corresponding to a given couple of 
parameters simply a model. 

As the fitted models are constructed such that their 
projected kinematics match these of the input models, we 
can talk about the projected kinematics of a model. The 
same applies to the light profile or quantities derived from 
it, such as the observed luminosity L b s (see section 2.2). On 
the contrary, when we describe spatial kinematic quantities 
(such as the anisotropy (3(r)) or integrals thereof (such as 
the mass-to- light ratio T), we need to distinguish between 
the ones corresponding to the input and fitted models, which 
a priori have no reason to be equal. With an apparent quan- 
tity we mean a quantity that corresponds to the fitted model, 
e.g. the quantity that results from the modelling procedure. 
With an intrinsic quantity we mean the quantity that corre- 
sponds to the input model, and hence is always independent 
of the optical depth. For example, all the models have the 
same intrinsic mass Mo, whereas the appar ent mass of the 
models is determined as outlined in section 3.3, and will be 
different for each model. Obviously, for optical depth r = 
the apparent and intrinsic values are equal. 

In this section we will compare the apparent and the 
intrinsic kinematic quantities of our models, as a function of 
the parameters r and q. 



4.1 The mass and the mass-to-light ratio 

In figure ^a we show the apparent dynamical mass of the 
models as a function of the optical depth. The global effect 
of the dust is clear : the mass decreases nearly linearly with 
the optical depth. The slope of this correlation however is 
strongly dependent on the orbital structure. For the very 



© 0000 RAS, MNRAS 000, 000-000 



6 Baes, Dejonghe & De Rijcke 



100 



10 r 



1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 



\ 

1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 



100 



- 10 r 



1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 



■ 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 



100 



- 10 r 



1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 



■ 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 



100 



1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 

p •••• 



1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 



2 4 6 8 10 2 4 6 8 10 02468 10 02468 10 

x (kpc) x (kpc) x (kpc) x (kpc) 



200 - 



in 

\ 

J 150 E 

a. 
b 



100 P. . . . . ,-3 ioo E- 

tXllluxuxluxllLlllI 



- 200 - 



,- 150 r 



-I 1 I I I 


|III|III|L 




1 1 1 1 1 1 | 
• 






1 1 1 1 1 




-1 1 1 1 1 1 1 1 1 1 1 1 1 


1 1 1 1 1 1 1 1 1 r 



- 200 - 



i 150 '- 



- 100 - 




- 200 - 



i 150 r 



- 100 - 




02468 10 02468 10 02468 10 02468 10 

x (kpc) x (kpc) x (kpc) x (kpc) 



-i i I i i i i I i i i i I i i. 




-500 , 500 
v (km/s) 



-500 , 500 
v (km/s) 



-500 t 500 
v (km/s) 



-500 t 500 
v (km/s) 



0.3 



o 
a 
a. 



o 0.2 - 



Q o.ih 

> 

o 
-J o I 



1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 




0.3 



- 0.1 - 



A 
/ \ 



0.3 



1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 
-I 0.2 h -I 0.2 - 



- 0.1 - 



1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 
: ys - 




0.3 



1 1 1 1 1 1 1 1 1 1 

1/ 


1 1 1 1 1 1 1 1 1 1 

\j 










- 





-400 400 -400 400 -400 400 -400 400 

v (km/s) v (km/s) v (km/s) v (km/s) 




- 250 - 



- 200 - 



- 150 - 




- 250 - 



- 200 



- 150 - 




02468 10 02468 10 02468 10 02468 10 

x (kpc) x (kpc) x (kpc) x (kpc) 



Figure 3. Results of the fitting procedure for the very tangential (left), tangential (middle left), isotropic (middle right) and radial 
(right) models, with optical depth r = 2. Shown are the projected light density profile t v (x) (in arbitrary units), the projected dispersion 
profile (J p (x) and the LOSVDs for x = Okpc and x = 10 kpc. The dots are the data points, the solid lines represent the fit. The bottom 
row shows the projected fourth moment £ p (x) of the LOSVDs, which is not included in the fit and can be used as a check on the results. 
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Figure 4. Three plots, showing the apparent mass M, the apparent mass-to-light ratio T and the apparent mean anisotropy B of our 
models, as a function of the optical depth t. They are shown for different orbital structures, ranging from very tangential (black) to 
radial (light grey). The dots represent the results of our fits, the solid lines are least-square fits to these points, either linear (for the 
mass and the mass-to-light ratio) or quadratic (for the mean anisotropy). The corresponding coefficients are tabulated in table bL 



Table 2. The parameters oaf, ax, ag and feg from the least- 
squares fits to AM, AT and A£> (see text). For hm and ax we 
give both the absolute and relative values. 



(10 9 M ) 

-6 0.57 (1.1%) 

-2 1.46 (3.0%) 

2.68 (5.4%) 

1 2.78 (5.6%) 



(T Q ) 

1.02 (25%) 
0.90 (23%) 
0.74 (19%) 
0.72 (18%) 



0.65 
0.24 
0.12 
0.02 



6 B 



-0.09 
-0.03 
-0.02 
-0.00 



tangential model the apparent mass is quite insensitive to 
the presence of dust, and the impact of the dust extinction 
becomes gradually stronger as we move to the tangential, 
the isotropic and the radial model. We fitted straight lines 
AM = aM t through the data to obtain characteristic values 
for the mass decrease in function of the optical depth. These 
values are given in the second column of table ^. 

The apparent luminosity of our models is calculated by 
integrating the apparent light density l(r) over space, or - 
since the fitted models contain no dust - by integrating l p (x) 
over the plane of the sky. It thus equals the observed lumi- 
nosity of the dusty Plummer models, which are tabulated in 
the third column of table 

Combining these with the apparent masses we can cal- 
culate the apparent mass-to-light ratio T of our models. The 
apparent luminosity decreases stronger than the apparent 
mass with increasing optical depth, such that T is an in- 
creasing function of r. The dependence of T on the orbital 
structure is only determined by the apparent mass, as the 
apparent luminosities are independent of the orbital mode. 
As a consequence, the dependence on q is now opposite: 
the mass-to-light ratio is most dramatically affected for tan- 
gential models, and the effects are smaller for radial and 
isotropic ones. The results are shown in figure ^b. Again, 
the dependence on the optical depth is nearly linear and 
straight lines AT = ar r are fitted through the data points 



to obtain characteristic values. These are tabulated in the 
third column of table 



4.2 The distribution function 

The eight panels in figure [E| represent the isoprobability con- 
tours of eight DFs, corresponding to the four t — 2 models. 
In the upper row we plot the intrinsic DFs, whereas the 
lower panels represent the apparent DFs. From left to right 
we have, as in the previous plots, the very tangential, tan- 
gential, isotropic and radial models. 

The contour plots are shown in turning point space, 
such that the DFs can easily interpreted in terms of or- 
bits. Let us first concentrate on the four upper panels. In 
the innermost regions of the galaxies, the shape of the DF is 
comparable - indeed all Plummer models are fairly isotropic 
in their centers. From a few kpc on however, we can clearly 
see how the isoprobability contours reflect the orbital struc- 
ture of the model they represent. Tangential models pre- 
fer nearly circular orbits, with a small difference between 
apocentre and pericentre, and their contours will tend to lie 
alongside the diagonal axis. On the other hand, radial mod- 
els prefer elongated orbits, with a large difference between 
apocentre and pericentre, such that their contours will tend 
to be more vertical. The slope of the isoprobability contours 
is thus indicative for the orbital structure of the model. 

Let us now compare the intrinsic and apparent DFs. 
Concerning the central regions we see that the same struc- 
ture is preserved for all orbital modes. Outside this region 
however, there are differences, most clearly visible for the 
very tangential and tangential models : the contours are ly- 
ing somewhat more horizontally, indicating that elongated 
orbits are relatively more favored. Dust obscuration thus 
seems to make these galaxies appear less tangential outside 
the innermost regions. Whether a similar trend accounts for 
the isotropic and radial models too, is less obvious from the 
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Figure 5. Isoprobability contour plots in turning point space for the r = 2 models. The upper and lower panels represent the intrinsic 
and the apparent DFs respectively. From left to right we have the very tangential, tangential, isotropic and radial models, as in figure pi 



DF contour plots. Therefore we will study the anisotropy of 
our models in detail. 



4.3 The anisotropy 

The intrinsic anisotropics of our models, as defined by equa- 
tion @, can be written as 



2 r 2 + c 2 



(18) 



where c represents the (true) core radius. All models are 
hence intrinsically isotropic in the central regions, at least 
if isotropy is defined only from the second order moments, 
and show their true orbital behavior at larger radii. 

The dependence of the apparent anisotropy on the op- 
tical depth is shown, for the different orbital modes, in fig- 
ure [| Shown are apparent (solid lines) and the intrinsic 
(dotted lines) anisotropies of the r = 1 and r = 3 models as 
a function of the spatial radius. For the very tangential and 
tangential models, we see that the effect, as determined from 
the DF plots, is confirmed: as for the intrinsic orbital struc- 
ture, all models are isotropic in their inner regions, and they 
are less tangential at larger radii. Rather logically, this effect 
increases with increasing optical depth. Looking at the third 
and fourth panel we see that also the apparent anisotropy of 
the isotropic and radial models increases outside a few kpc. 
All models thus seem to be subject to a "radialization" , i.e. 
dust obscuration tends to make galaxies appear more radi- 
ally anisotropic outside the central few kpc, even if they are 
already intrinsically radial. 

In order to quantify the strength of the radialization 
in function of the input parameters q and r, it is useful to 



consider one single anisotropy parameter. We define a mean 
anisotropy B as 



B = 



J p(r) p(r) r 2 dr 



' dr 



(19) 



with the integral covering the region of our fits (between 
and 10 kpc). The intrinsic mean anisotropy for our models is 
off course proportional to the parameter q; substituting ( |l8| ) 
and (0) in M) we find B = ^ q. 

In figureTjc we plot the apparent mean anisotropy B 
for our models as a function of the optical depth. One can 
clearly see that B increases for increasing optical depth, and 
that the radialization is more dramatic the more tangential 
the input model. The curves in the figure are quadratic fits 
AB — asT + brsT 2 to the data points, and the coefficients 
are tabulated in the last two columns of table 0. 



5 DEPENDENCE ON THE DUST MODEL 

The calculations in the previous chapters are based on a dust 
model whose spatial dependence is given by equation (^). 
But, as discussed in the introduction, very little is known 
about the spatial distribution of the dust in elliptical galax- 
ies. In this chapter we will investigate whether the results so 
far obtained change dramatically if the relative distribution 
of dust and stars varies. Therefore we consider, as in paper I, 
a more general opacity function, 



n(r) — —= r 



i r(f) 



' r (sESJ c{ + c 2 



(20) 



which also satisfies the normalization condition (g), and 
which reduces to (El) if a = 3. The extra parameter in this 
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Figure 6. The orbital structure of input and output models, as characterized by the anisotropy From left to right we have a panel 
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Figure 7. The opacity function n(r), corresponding to equa- 
tion j2o| ) for different dust exponents a. It is shown for r = 1, 
for other values of r it scales linearly. For large values of a, the 
dust is concentrated in the central regions, for large values the 
distribution is extended. 



family of opacity functions, the dust exponent a, sets the 
spatial distribution of the dust. In figure M we plot the opac- 
ity function for different values of a. 

Small values of a correspond to spatially extended dust 
distributions. The range of a is restricted to a > 1. When 
a approaches this value, the dust is more or less equally 
distributed along the line-of-sight, and as c <C D, relatively 
very little dust resides in the central regions of the galaxy. 
In the limit a — > 1 the opacity function is such that the dust 
effectively forms an obscuring medium between the galaxy 
and the observer (see paper I), analogous to extinction of 
starlight due to interstellar dust in the Galaxy. This geom- 
etry is generally known as the overlying screen approxima- 
tion, and it is the geometry which has, for a fixed optical 
depth, the largest impact on the projection of starlight. For 
many years the extinction in spiral galaxies was described 
using the cosecans law (Holmberg 1975), which implicitly as- 
sumes this geometrical distribution of stars and dust. Nowa- 
days however, the extinction in these systems has been de- 
scribed using more detailed dust-stars geometries, and the 



Table 3. Same as tabic |l| but now for a fixed optical depth r = 2 
and as a function of the dust exponent a. The last row, which is 
labeled ND, shows the same quantities for the dust-free Plummer 
models. 



Q 


'-obs 


^obs 


A 


?p,0 


a P,0 




(kpc) 


(10 9 L ) 


(mag) 


(L /pc 2 ) 


(km/s) 


1.0 


5.00 


4.60 


1.09 


58.6 


178.0 


1.5 


5.58 


5.77 


0.84 


59.7 


177.6 


2.0 


5.91 


6.75 


0.67 


61.5 


177.1 


3.0 


6.35 


8.11 


0.47 


64.6 


176.5 


5.0 


6.69 


9.55 


0.29 


68.8 


176.1 


ND 


5.00 


12.50 


0.00 


159.2 


178.0 



overlying screen model is generally considered to be unsat- 
isfying (Bruzual, Magris & Calvet 1988; Disney, Davies & 
Phillips 1989; Witt, Thronson & Capuano 1992). 

On the other hand, larger values of a correspond to 
centrally concentrated dust. For a = 5 dust and stars have 
the same geometry, and if a becomes very large the extinc- 
tion is confined to the central regions of the galaxy only. 
Silva & Wise (1996) investigated the effects of centrally con- 
centrated dust distributions on the photometry of elliptical 
galaxies. They found that, for models where the stars and 
dust have the same spatial distribution or where the dust 
is more concentrated than the stars, steep color gradients 
would be implied in the core, even for small optical depths. 
However, Crane et al. (1993) and Carollo et al. (1997) im- 
aged the cores of a set of nearby elliptical galaxies using 
HST, and both of them found relatively small color gradi- 
ents, and hence no direct sign for the presence of centrally 
concentrated diffuse dust distributions. 

We consider a = 1.5, a = 2 and a — 5 (besides a — 3) 
and create new dusty Plummer models for each of these ex- 
ponents and for the four orbital modes, where we fix the 
optical depth at the median value r = 2. Data sets are cre- 
ated and these are modelled exactly as before. In particular, 
the observed light profile can still in a satisfying way be ap- 
proximated by a Plummer potential (|^), for all values of a 
under consideration. The observed core radii c b s , as well 
as some other parameters of the input models, are listed in 
table | as a function of the dust exponent. In figure ^ we 
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Figure 8. The apparent mass M, the apparent mass-to-light ratio T and the apparent mean anisotropy B, as a function of the dust 
exponent a, with r = 2 fixed. They are shown for the different orbital structures, ranging from very tangential (black) to radial (light 
grey). The dots represent the calculated values, the solid lines are just lines to guide the eye, in contrary with figure W. The asterisks 
represent the intrinsic values of the represented quantities. 



plot the results of our modelling: the three plots show the 
apparent dynamical mass M, the apparent mass-to-light ra- 
tio T and the apparent mean anisotropy B as a function 
of the dust exponent a, for the four different values of q. 
The intrinsic values of these quantities are indicated by as- 
terisks, in contrast with figure ^j, where these correspond to 
the r = case. 

Qualitatively, the effects are the same, independent of 
the value of the dust exponent: the apparent dynamical mass 
decreases, the orbital structure seems more radial and the 
apparent mass-to-light ratio increases. Moreover, the depen- 
dence on the orbital structure of these effect is independent 
of the dust exponent a: the radial and isotropic models tend 
to lose more mass and keep their orbital structure, whereas 
for the tangential and very tangential models it is vice versa. 

Quantitatively, there is a dependence on the dust geom- 
etry. The apparent dynamical mass and the apparent orbital 
structure are unaffected for both very large or very small 
values of a. The effects appear to be strongest for a around 
3. The mass-to-light ratio however, has totally different de- 
pendence on a, it is strongly affected (a factor 2 to 3) for 
extended dust distributions, and the effects decreases grad- 
ually when the dust becomes centrally concentrated. In fact, 
T is largely determined by the apparent luminosity, which is 
a strong function of the dust geometry: the extinction is far 
more effective for extended than for condensed distributions 
(paper I, WS96). 

This behavior can be illustrated when we consider the 
limits of very low and very high dust exponents. For an 
extremely extended dust distribution (a — * 1) on the one 
hand, the dust effectively forms an absorbing screen with 
optical depth J between the galaxy and the observer. Rel- 
ative to the dust-free case, the light profile then decreases 
with a factor exp f £), independent of the TOS, whereas the 
morphology of £ P (x) and the shape of the LOSVDs are not 
affected (paper I). Hence, a model with exactly the same 
potential, dynamical mass and orbital structure, but with a 



mass-to-light ratio which is a factor exp higher, will fit 
these data exactly. For an extremely centrally concentrated 
dust distribution (a ^> 5) on the other hand, the dust effects 
will be visible only in the very innermost regions. As a result, 
neither the projected light density, nor the LOSVDs will be 
severely affected, and hence the data will be nearly the same 
as in the dust-free model. Hence all intrinsic quantities are 
recovered. 

The question we want to answer in this chapter was 
whether the results we obtained using the a = 3 model 
change dramatically in function of the dust geometry. As 
we argued, neither centrally concentrated nor very extended 
dust distributions seem very probable. Therefore we con- 
sider the range 2 55, a 5, where the stars are somewhat 
more concentrated than the dust, as representative for re- 
alistic geometries. Although M, B and T do vary with a 
in this range, the effects have qualitatively the same behav- 
ior (decreasing apparent dynamical mass, increasing appar- 
ent mean anisotropy), and quantitatively the same order of 
magnitude. Hence we can conclude that our results of the 
a — 3 case, as summarized in table can be considered as 
representative. 



6 DISCUSSION 

In this paper we investigated which errors can be made by 
not taking dust into account in dynamical modelling proce- 
dures. Therefore we created a set of galaxy models consist- 
ing of a dust and a stellar component. We calculated the 
projected kinematics, taking dust into account, using the 
method outlined in paper I. These data sets are then mod- 
elled as if no dust were present, and the apparent dynamical 
properties of these models are calculated and compared to 
the intrinsic ones, as a function of the orbital structure of 
the input model, the optical depth of the dust and the dust 
geometry. 

We find that (1) the dynamical mass of the galaxy tends 
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to become smaller, (2) the orbital structure seems to be radi- 
alized. For a fixed optical depth and dust geometry, the rela- 
tive strength of these effects depends on the orbital structure 
of the input model. For radial and isotropic models the ap- 
parent mass decreases significantly, with a typical amount 
of 5 per cent per optical depth unit, whereas their orbital 
structure is hardly affected. For tangential models on the 
other hand, the dynamical mass is less sensitive to the pres- 
ence of dust (about 3 per cent per optical depth unit for 
the tangential model and less for the very tangential model), 
whereas now the radialization is considerably stronger. Both 
effects are apparently coupled, in a way that reminds of the 
mass-anisotropy degeneration in spherical systems. 

The effects are dependent on the shape of the dust dis- 
tribution, but not very critically: for the dust exponents in 
the range 2 £ a £ 5 the effects are very comparable. On the 
one hand this is fortunate, since it means that our calcula- 
tions are more or less model-independent, and can be appli- 
cated for a wide range of dust geometries. On the other hand 
however, this means that dynamical analyses will not be able 
to discriminate convincingly between different values for the 
dust exponent. Analyses as ours can thus hardly be used to 
constrain the distribution (and the origin) of the smooth 
dust component in ellipticals. Further infrared and submil- 
limeter observations, in particular ISO data, are necessary to 
solve this problem. Preliminary results include the detection 
of an extended, very cold dust component in the dwarf ellip- 
tical NGC205 (Haas 1998), and of warm dust in the central 
regions of the Seyfert I SO galaxy NGC 3998 (Knapp et al. 
1996). A larger database of ISO imaging of early- type galax- 
ies at both mid-infrared an far-infrared wavelengths would 
improve our knowledge significantly. 

The combined results of paper I and this paper may at 
first glance seem quite contradictory. In paper I we found 
that the observed kinematics of elliptical galaxies are not 
severely affected by dust obscuration. Whence it seemed ob- 
vious that modest amounts of dust do not imply large un- 
certainties on dynamical mass determinations or estimates 
of the anisotropy of these systems. In this detailed study 
however we find that dust does have an important effect on 
the determination of the dynamical structure, in particu- 
lar the dynamical mass and the anisotropy. The answer to 
this apparent discrepancy lies in the fact that the potential 
plays am important role in the determination of the inter- 
nal structure of galaxies. Dejonghe & Merritt (1992) show 
that, in case of a spherical two-integral system, the knowl- 
edge of the potential and the entire set of LOSVDs suffice to 
determine the DF uniquely. The potential itself can be con- 
siderably constrained by the LOSVDs, but is not uniquely 
determined, such that a set of potentials will usually yield 
acceptable models for a data set. Often one chooses that po- 
tential that is derived from the observed light profile, if this 
one is one of the possible choices (e.g. if dark matter is not 
assumed to play a major role). But if dust is present, the 
light profile will be severely affected, even by small amounts 
of dust, such that the matching potential will not be the 
correct one. And as diffuse dust is assumed to be present in 
a major fraction of the early-type galaxies, we argue that it 
is important to at least be aware of its effects, which may 
not be as trivial as one might imagine. As WS96 stressed, 
all broadband observations of elliptical galaxies may be af- 
fected by dust, and hence dust should be seriously taken into 



account in their interpretation. We now can add that dust 
does also play a role in dynamical analyses, and hence that 
it should also here be taken into account, in a non-trivial 
way. 

We close by giving a simple example to illustrate this 
point. A simple way to estimate the mass of a gravitating 
system is its virial mass. For example, Tonry & Davis (1981) 
estimate the masses for a set of 373 elliptical galaxies using 
a relation where mass is proportional to the effective ra- 
dius and the square of the central dispersion. Although it is 
nowadays possible to obtain much better mass estimates for 
nearby galaxies, the virial mass estimate is still one of the 
only tools to constrain the mass of galaxies and clusters at 
intermediate or high redshifts (Carlberg et al. 1996, Carl- 
berg, Yee & Ellingson 1997, Tran et al. 1999). The question 
is now how to correct these mass estimates for the presence 
of dustQ. A straightforward way is to estimate amount of 
dust using IRAS or ISO data, and calculate the effects on 
dispersion and scale length. Dispersions are only slightly af- 
fected by dust absorption (paper I), whereas scale lengths 
as the effective radius or the core radius can increase sub- 
stantially, as dust primarily removes light from the centre 
of the system. Hence we find that the apparent mass of the 
galaxy would increase as a function of r, while we find, us- 
ing detailed kinematic modelling, that the apparent mass 
decreases with increasing optical depth. Moreover, this cor- 
rection will be independent of the orbital structure of the 
model. This again illustrates the fact that dust effects are 
non-trivial and should be fully taken into account. 
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